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Abstract 

We construct an explicit, TeV-scale model of decaying dark matter in which the approximate 
stability of the dark matter candidate is a consequence of a global symmetry that is broken only 
by instanton-induced operators generated by a non-Abelian dark gauge group. The dominant dark 
matter decay channels are to standard model leptons. Annihilation of the dark matter to standard 
model states occurs primarily through the Higgs portal. We show that the mass and lifetime of 
the dark matter candidate in this model can be chosen to be consistent with the values favored by 
fits to data from the PAMELA and Fermi LAT experiments. 
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I. INTRODUCTION 



Evidence has been accumulating for an electron and positron excess in cosmic rays com- 
pared with expectations from known galactic sources. Fermi LAT [l| and H.E.S.S. [2] have 
measured an excess in the flux of electrons and positrons up to a TeV or more. The PAMELA 
satellite is sensitive to electrons and positrons up to a few hundred GeV in energy, and is able 
to distinguish positrons from electrons and charged hadrons. PAMELA detects an upturn 
in the fraction of positron events beginning around 7 GeV [3|. This is in contrast to the 
expected decline in the positron fraction from secondary production mechanisms. Curiously, 
no corresponding excess of protons or antiprotons has been detected ^]. 

Although conventional astrophysical sources may ultimately prove the explanation of the 
anomalous cosmic ray data 5|, an intriguing possibility is that dark matter annihilation or 
decay provides the source of the excess leptons. If dark matter annihilation is responsible for 
the excess leptons, then the annihilation cross section typically requires a large boost factor 
~ 100 — 1000 to produce the observed signal [6j. Possible sources of the boost factor include 
Sommerfeld enchancement from additional attractive interactions in the dark sector {7], 



WIMP capture 



9] or Breit-Wigner resonant enhancement [10Ml2] . 
Alternatively, decaying dark matter can provide an explanation of the cosmic ray data 
if the dark matter decay channels favor leptonic over hadronic final states [13]. A typical 
scenario of this type that is consistent with PAMELA and Fermi LAT data includes dark 
matter with a mass of a few TeV that decays to leptons, with an anomalously long lifetime of 
~ 10 26 seconds From a model-building perspective, an intriguing issue is the origin of 

this long lifetime, and whether it can be explained with a minimum of theoretical contrivance. 
With this goal in mind, we present a new model of TeV-scale dark matter, one in which an 
anomalous global symmetry prevents dark matter decays except through instantons of a non- 
Abelian gauge field in the dark sector. Instanton-induced decays naturally produce the long 
required lifetime. Small mixings between standard model leptons and dark fermions gives 
rise to the leptonic final states observed in the cosmic ray data. Dark matter annihilation 
through the Higgs portal allows for the appropriate dark matter relic abundance, with dark 
matter masses consistent with the range preferred by PAMELA and Fermi-LAT data. 

Superheavy dark matter decays through instantons have been considered before as a 
possible explanation for ultra-high energy cosmic ray signals, but those scenarios assumed 
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FIG. 1: Dark matter decay vertex. The circle represents the instanton- induced interaction, while 
X's represent mass mixing between the \ fields and standard model leptons. Note that e and v 
represent leptons of any generation. 



superheavy dark matter with a mass of 10 13 GeV or higher 16( which cannot simultane- 
ously explain the lower energy electron and positron flux being considered here. Models of 
anomaly-induced dark matter decays without a dark gauge sector can also be constructed. 
For example, a supersymmetric extension of the radiative seesaw model of neutrino masses 
can explain the PAMELA data through dark matter decays via an anomalous discrete sym- 
metry 17J . The TeV-scale model we present, which is based on the smallest, continuous 
non-Abelian dark gauge group and smallest set of exotic particles necessary to implement 
our idea, suggests a prototypical set of new particles and interactions that could perhaps be 
probed at the LHC. 

In Section [TTlwe present the model and describe the leptonic decay mode via instantons. In 
Section HTT1 we consider dark matter annihilation channels and demonstrate that annihilation 
through the Higgs portal can lead to the measured dark matter relic density. In Section IIVI 
we consider dark matter interactions with nuclei and find that our model is safely below 
current direct detection bounds. We conclude in Section [V] 
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#g (2,0)o 77 (l,l/6)o 

TABLE I: Particles charged under the dark gauge groups. The SU(2)£>xU(l)£) charge assignments 
are indicated in parentheses; the subscripts +, — and represent the standard model hypercharges 
+1, —1 and 0, respectively. Note that the if) and x states are fermions, while the Hp and r\ are 
complex scalar s. 

II. THE MODEL 

The gauge group of the dark sector is SU(2)cxU(1)d. The matter content consists of 
four sets of left-handed SU(2)d doublets and right-handed singlets: 

Vfc = j i>uR, i^m ; X? = ^ ^ j xl, xf R (< = l • • • 3) (2.1) 

We include an SU(2) D doublet and singlet Higgs field, H D and rj, respectively, that are 
responsible for completely breaking the dark gauge group. In addition, the Higgs field H D 
is responsible for giving Dirac masses to the if) and \ fields. The model is constructed 
so that if) number corresponds to an anomalous global symmetry that is violated by the 
X^X vertex generated via SU(2)£> instantons, as indicated in Fig. [TJ The x fields 
are assigned hypercharges so that they mix with standard model leptons, leading to the 
decay if) — > i + i~v. The required lifetime (~ 10 26 s) and the appropriate dark matter relic 
density [Vtj^h 2 ~ 0.1) constrain the free parameters of the model. 

The charge assignments for these fields are summarized in Table [H Let us first discuss 
the consistency of the charge assignments. Cancellation of the SU(2)|, U(l) anomalies 
requires that the sum of the U(l) charges over all the dark doublet fermion fields must 
vanish. As one can see from Table [H this is clearly the case for the U(1)d and U(l)y 
charges of the left-handed doublet if) and x fields. Since SU(2) is an anomaly free group 
and has traceless generators, all other SU(2) r> anomalies vanish trivially. Now consider the 
U(l)^,U(l) y anomalies (where p and q are non-negative integers satisfying p + q = 3). For 
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each field in Table |T] with a given U(l)oxU(l)y charge assignment, one notes that there is 
another with the same charge assignment but opposite chirality. As far as the Abelian groups 
are concerned, the theory is vector-like and the corresponding anomalies vanish. Finally, we 
note that the theory has precisely four SU(2)o doublets and is free of a Witten anomaly. 

The gauge symmetries of the model lead to a global U(l)^, symmetry that prevents the 
decay of the lightest if) mass eigenstate at any order in perturbation theory. To confirm this 
statement, we need to show that all renormalizable interactions that violate this symmetry 
are forbidden by the dark-sector gauge symmetry. The possible problematic interactions 
that could violate this global symmetry fall into the following categories: 

1. Terms involving if> c ifj. Here the superscript indicates charge conjugation, ijj c = i^^ 2 if) T . 
This combination has U(l)^ charge +2. However, it also has U(1)d charge —1. Since we 
have no Higgs field with the U(1)d charge ±1, there are no renormalizable interactions that 
violate if> number by two units. 

2. Terms involving a x fermion and if) or ip c . Such terms violate if) number by ±1 unit. 
However, the possible bilinears involving if) and any \ have U(1)d charges ±1/3 or ±2/3. 
Again, we have no Higgs field with the necessary U(l)o charge to form a renormalizable 
gauge invariant term of this type. 

3. Terms involving a standard model fermion and if) or ip c . Such an interaction would 
violate if) number by ±1, but would have U(1)d charge ±1/2. Again, we have no Higgs 
fields with charge ±1/2 that would allow the construction of a renormalizable invariant. 

Since the renormalizable interactions of the theory have an unbroken U(l)^, symmetry, no 
perturbative process involving these interactions will violate the global symmetry. However, 
since the SU(2)f, U(l)^ anomaly is non-zero, non-perturbative interactions due to instantons 
will generate operators that violate the U(l)^ symmetry. 

Instantons are gauge field configurations which stationarize the Euclidean action but have 
a nontrivial winding number around the three-sphere at infinity. Following 't Hooft 18, 



191 ] . if there are Nf Dirac pairs of chiral fermions which transform in the fundamental 



representation of a gauge group, then due to the chiral anomaly a one-instanton configuration 
violates the axial U(l)^ charge by 2Nf units. The non- Abelian, SU(Nf)xSU(Nf) chiral 
symmetry is non-anomalous, so the instanton process must involve the 2Nf chiral fermions 
in a symmetric fashion. Fig. [1] shows the effective ipX^X^X^ interaction induced by the 
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instanton configuration in our model. 1 Given the hypercharge assignments of the x fields, 
these states have electric charges +1, and —1, the same as standard model leptons, of 
any generation. After the dark and standard model gauge symmetries are spontaneously 
broken, there is no symmetry which prevents the x states and the standard model leptons 
from mixing. By including a single vector-like lepton pair, we now show that mixing leading 
to the decay ip ~ >* can arise via purely renormalizable interactions. 

We introduce a vector-like lepton pair, El, Er, with mass Me and the same quantum 
numbers as a right-handed electron; in the notation of Table [TJ 

El ~ E/r ~ (1, 0)_ . (2.2) 

In addition, we assume in this model that standard model neutrinos have purely Dirac 
masses. If the Higgs vacuum expectation values (vevs) are smaller than the masses of the 
heavy states, then the mixing to standard model leptons shown in Fig.[T]can be estimated via 
the diagram in Fig. |2j Otherwise, one has to diagonalize the appropriate fermion mass ma- 
trices. We discuss the exact diagonalization in an appendix for the reader who is interested 
in the details. Here, the diagrammatic approach is sufficient to establish that the mixing is 
present, and is no larger than order (r))/M x , (r])/M x , and (rj) (H) / \M X M E ) , where H is the 
standard model Higgs, for the x^l ~ e Ri Xl? ~~ v r anc ^ Xl ~ e L mixing angles, respectively. 
We take each mixing angle to be 0.01 in the estimates that follow, and demonstrate in the 
appendix how this choice can be easily obtained. Further, we assume that decays to the 
heavy eigenstates are not kinematically allowed, as is also illustrated in the appendix. Due 
to the mixing, the particles decay quickly to standard model particles via couplings 
to the Higgs bosons and standard model electroweak gauge bosons. The heavier ip mass 
eigenstate decays to lighter states via SU(2)d gauge-boson-exchange interactions. 
The instanton-induced vertex in Fig. [1] follows from an interaction of the form 

C ( 8tt 2 \ /m^\ 35/6 1 
6g D V 9d J \ vdJ v z d 

6S>2)(^xi 3)7 ) - (^2)(^!?xl 3)7 )l + h.c. , (2.3) 



where a, (3, 7 and a are SU(2)£> indices |l9|, |20j. The dimensionless coefficient C can 



be computed using the results in Ref. 19| and one finds C ~ 7 x 10 8 . The operators in 



1 In this model, Planck-suppressed operators of this form, if they are present, are negligible compared to 
the instanton-induced effects. 
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FIG. 2: Diagrammatic interpretation of mixing from \ states to standard model fermions, corre- 
sponding to the right-hand-side of Fig. [TJ Here E represents the vector-like lepton described in the 
text, and H is the standard model Higgs. 

Eq. (I2.3P lead, via mixing, to operators of the form vr^l^r^l and (^b^lPr^-l- Assuming 
that the product of mixing angles is ~ 10~ 6 , as discussed earlier, one may estimate the decay 
width: 



For example, for = 3.5 TeV and vd = 4 TeV, one obtains a dark matter lifetime of 10 26 s 
for 



where go is defined in dimensional regularization and renormalized at the scale [19J. 
For similar parameter choices, one can slightly adjust to maintain the desired life- 
time. As mentioned earlier, dimension-six Planck-suppressed operators are much smaller 
than the operators in Eq. (I2.3p . Sphaleron- induced interactions are suppressed by ~ 
exp[— 4ttv d / (goT)} ~ exp(— 44 TeV/T), and become negligible well before the tempera- 
ture at which dark matter freeze out occurs. 

Finally, let us consider whether the choice = 4 TeV conflicts with other meaningful 
constraints on the heavy particle content of the model. In short, a spectrum of ~ 4 TeV 
X and E fermions with order 0.01 mixing angles with standard model leptons presents no 
phenomenological problems. These states are above all direct detection bounds; they are 




(2.4) 



9d ~ 1-15, 



(2.5) 



7 



vector-like under the standard model gauge group so that the S parameter is small; they mix 
weakly enough with standard model leptons so that other precision observables are negligibly 
affected. On this last point, we note that the correction to the muon and Z-boson decay 
widths due to the fermion mixing is a factor of 10~ 8 smaller than the widths predicted in 
the standard model, which is within the current experimental uncertainties. The dark sector 
gauge bosons are also phenomenologically safe. They do not have couplings that distinguish 
standard model lepton flavor (since they do not couple directly to standard model leptons) 
so that tree- level lepton-flavor violating processes are absent. The effective four-standard- 
model-fermion operators that are induced by dark gauge boson exchanges are suppressed by 
~ (0.01) 4 /t )p ~ 1/(40,000 TeV) 2 , which is consistent with the existing contact interaction 



bounds (25] 



We now turn to the question of whether the model provides for the appropriate dark 
matter relic density. 

III. RELIC DENSITY 

For the regions of model parameter space considered in this section, dark matter annihi- 
lations to standard model particles proceed via mixing between the dark and ordinary Higgs 



bosons, often described as the Higgs portal 



21] . We take into account mixing between the 



doublet Higgs fields, Hp and H, in our discussion below. This is consistent with a simplify- 
ing assumption that the r\ Higgs does not mix with the others in the scalar potential. Such 
an assumption is adequate for our purposes since we aim only to show that some parameter 
region exists in which the correct dark matter relic density is obtained. Consideration of a 
more general potential would likely provide additional solutions in a much larger parameter 
space, but would not alter the conclusion that the desired relic density can be achieved. 

In this section, if) will refer to the dark matter mass eigenstate, i.e., the lightest mass 
eigenstate of the if) u -if>d mass matrix, which we take as diagonal, for convenience. The 
potential for the doublet fields has the form: 

V = -^H^H + X(H^H) 2 - filH^Hu + \ D {H ] D H D f + X mix (H^ H)(H^ D H D ). (3.1) 



S 



In unitary gauge, H and Hp are given by 



V2 





v + h 



H D = ^ 



V2 





vd + h D 



(3.2) 



where v and vd are the H and Hd vevs, respectively. At the extrema of this potential, 

1 



V (-/J 2 + X V 2 + -Xmix V 2 D ) 
V D {~iAd + V 2 D + ~X m ix V 2 ) 











(3.3) 



The h-ho mass matrix follows from Eq. (13. ip . 



M 2 H 



2Xv 2 X 



vv D 2X D v 2 




(3.4) 



Diagonalizing the mass matrix, one finds the mass eigenvalues 



m 



1,2 



(X D v 2 D + Xv 2 )t (Xdv 2 d - Xv 2 )^l + ^\ 



where 



V 



XmixV v D 



X D v 2 D - Xv 2 ' 

The mass eigenstates hi and hi are related to h and /io by a mixing angle 

h\ = hcosO — hpsmO 
hi = h sin 9 + hp cos 9, 



(3.5) 



(3.6) 



where 



tan 2 9 = y . 



(3.7) 



(3.* 



Dark matter annihilations proceed via exchanges of the physical Higgs states hi and 
h 2 . We take into account the final states W + W~, ZZ, h-Jii and tt, where t represents 
the top quark. For the parameter choices considered later, final states involving hi will be 
subleading. The relevant annihilation cross sections are given by 



g 2 mL sin 2 9 cos 2 9 2 
&w+w- — — tt^, 2 — 2 s ' 



1 



1 



s — mj + imiTi s — m\ + zm 2 r 2 



xi 1 



4m 2 / in i 



(3.9) 



9 



&zz 



g 2 m^ sin 2 9 cos 2 9 2 
2bQ^m 2 v v 2 D 



m\ 




Am 2 , 



4m| 



1 - 



irriiTi 
Am 2 



s — m\ + im-iY-i 



(3.10) 



16^1, 



3#msin0 



5-112 cos 



s — m 2 + imiTi 



s — mo + im2^2 




Am\ 



Am 2 hi 



(3.11) 



3m 2 m 2 sin 2 cos 2 



Am 2 



1 



s — mj + imiTi s — m\ + ot^I^ 



4m 2 



( 3 - 12 ) 

s J \ s 

In Eqs. f )3.9p and ( I3.10p . g is the standard model SU(2) gauge coupling. In Eq. (13. 1 II) . gm 
and 5n2 represent the h\ and /i 2 ^i couplings, respectively: 



Qui = (Acos 3 + -A mia; cos0sm 2 i 



(Xd sin 3 + -A mi:E sin cos 2 0) t> D 



#112 = [3 A cos 2 6 1 sin — X mix (cos 2 sin — - sin 3 0)] t> 

+ [3 Xd sin 2 cos — X m i X (sin 2 6* cos — — cos 3 1 



(3.13) 



Finally, in all our annihilation cross sections, Fx (T 2 ) represents the decay width of the 
Higgs field hi (/i2)- The width Ti is comparable to that of a standard model Higgs boson 
and can be neglected without noticeably affecting our numerical results. However, since our 
eventual parameter choices will place the mass of the heavier Higgs field around 2m^, we 
must retain I^; the leading contributions to I^ come from the same final states relevant to 
the ip annihilation cross section: 



h 2 ->W+W- 



r 



2 3 

9 m% 
QATim^y 

2 1 
9 ™2 

l2Mm 2 w 

9n2 

327rm2 
3m2m 2 



sin 9 a 



sin 9\ 



4mfy 

2 

mi 



Am 2 z 

2~~ 

mi 



2 

mi 



+ 



1 - 



4m? 



12m^ 
mi 



I2m% 
mi 



Ami 



m% 



sin 2 9 1 



Am 2 
m 2 , 



3/2 



(3.14) 
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m^TeV) \/2Aw2(TeV) J 2\ D v 2 D (TeV) \ mix mi(GeV) m 2 (TeV) 



1.0 


0.19 


1.98 


0.21 


158 


1.98 


1.5 


0.22 


2.98 


0.28 


199 


2.98 


2.0 


0.26 


3.97 


0.39 


241 


3.97 


2.5 


0.27 


4.97 


0.42 


257 


4.97 


3.0 


0.29 


5.96 


0.52 


277 


5.96 


3.5 


0.31 


6.96 


0.57 


299 


6.96 


4.0 


0.35 


7.95 


0.70 


339 


7.95 



TABLE II: Examples of viable parameter sets for vp = 4 TeV. For each point listed, f^/i 2 ~ 0.1 
and the Higgs masses are consistent with the LEP bound. 

The evolution of the ip number density, n^, is governed by the Boltzmann equation 

^ + mt)n, = -{av)[n% - « Q ) 2 ], (3.15) 



where H(t) is the Hubble parameter and is the equilibrium number density. 



'he 



thermally-averaged annihilation cross section times relative velocity (cry) is given by [22] 

If 00 

where a tot is the total annihilation cross section, and the Ki are modified Bessel functions 
of order i. We evaluate the freeze-out condition [23] 

r _ n^ Q {av) 

mih) = ~WfV * 1 ' (3 - 17) 

to find the freeze-out temperature Tf, or equivalently Xf = m^/Tf. We assume the non- 
relativistic equilibrium number density 

3/2 



^ = 2 (l^) e_mWT ' (3 ' 18) 

and the Hubble parameter H = 1.66 g l J 2 T 2 /mpi, appropriate to a radiation-dominated 
universe. The symbol g* represents the number of relativistic degrees of freedom and rripi = 
1.22 x 10 19 GeV is the Planck mass. For the parameter choices in Tables HT1 and HTT| we find 
Xf ~ 27-28. We approximate the relic abundance using [2^] 



xo „i/a 



Y = Y s + V -k mpim ^ J Xf ,g 1^ {av)dx (3 ' 19) 
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m^TeV) \/2Aw2(TeV) J 2\ D v 2 D (TeV) \ mix mi(GeV) m 2 (TeV) 



1.0 


0.16 


1.98 


0.21 


121 


1.98 


1.5 


0.15 


2.98 


0.28 


118 


2.98 


2.0 


0.16 


3.97 


0.39 


127 


3.97 


2.5 


0.15 


4.97 


0.42 


124 


4.97 


3.0 


0.15 


5.96 


0.52 


122 


5.96 


3.5 


0.15 


6.96 


0.57 


127 


6.96 


4.0 


0.15 


7.95 


0.70 


122 


7.95 



TABLE III: Examples of viable parameter sets for vr> = 4 TeV, with mi below 130 GeV. For each 
point listed, Sl/j/i 2 ~ 0.1 and the Higgs masses are consistent with the LEP bound. 

where Y is the ratio of the number to entropy density and the subscript indicates the 
present time. The ratio of the dark matter relic density to the critical density p c is given by 
Qe> = Y^sqiti^/ p c , where s is the present entropy density, or equivalently 

Vt D h 2 « 2.8 x 10 8 GeV" 1 Y . (3.20) 

In our numerical analysis, we assume that the heavy states are sufficiently nondegenerate, 



24j. In Tables IT] and UTIl we 



so that we do not have to consider co-annihilation processes 
show representative points in the model's parameter space, spanning a range of ijj masses, 
in which we obtain the correct dark matter relic abundance, fljjh 2 w 0.1, and in which the 
masses m\ and m 2 are consistent with the LEP bound m 1>2 > 114.4 GeV |25|. It is common 
wisdom that weakly interacting dark matter candidates with masses of a few hundred GeV 
typically yield relic densities in the correct ballpark. We have assumed masses above 1 TeV 
since most fits to the positron excess in PAMELA and Fermi LAT indicate that a decaying 
dark matter candidate should have a mass in this range. One would therefore expect that 
Qoh 2 in our model should be larger than desirable. The reason this is not the case is that 
we have chosen parameters for which the heavier Higgs h 2 is within 1% of 2m^, leading 
to a resonant enhancement in the annihilation rate. While we would be happier without 
this tuning, it is no larger than tuning that exists in, for example, the Higgs sector of the 
minimal supersymmetric standard model. It is also worth pointing out that this tuning 
is related to the portal that connects the dark to standard model sectors of the theory 
and is not strictly tied to the mechanism that we have proposed for dark matter decay. 
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Other portals are possible. For example, one might study the limit of the model in which 
the U(1)_d gauge boson is lighter and kinetically mixes with hypercharge, a possibility that 
would lead to other annihilation channels. Finally, we point out that Tables [TT] and III II 
includes = 3.5 TeV, which naively corresponds to the value preferred by a fit to the 



PAME 



A and Fermi-LAT data, assuming a spin-1/2 dark matter candidate that decays to 



151 ]. However, other masses should not be discounted since astrophysical sources may 



also contribute to the observed positron excess 



IV. DIRECT DETECTION 



We now consider whether the parameter choices described in the previous section are 
consistent with the current bounds from direct detection experiments. The most relevant 
constraints come from experiments that search for spin- independent, elastic scattering of 
dark matter off target nuclei. The relevant low-energy effective interaction from t-channel 
exchanges of the Higgs mass eigenstates is given by 



it it 



where 



m q m^ sin 9 cos 9 
vv D 



qq 



1 



1 



VTLr 



(4.1) 



(4.2) 



This interaction is valid for momentum exchanges that are small compared to mi^, which 
is always the case given that typical dark matter velocities are non-relativistic. Following 
the approach of Ref. 



23], Eq. (14. ip leads to an effective interaction with nucleons 

£e// = fp $ippp + f n ip4> nn , 



(4.3) 



where f p and f n are related to a„ through the relation 



27] 



m 



p. it 



E 

q=u,d,s 



fp,n _ ^Tq °g 2 ( pjn ) 

~r ~- J' 



rn' 

q=c,b,t q 



(4.4) 



where (n\m q qq\n) = m n f^ q . Numerically, the f^^ 1 are given by 22] 

f Tu = 0.020 ± 0.004, f Td = 0.026 ± 0.005, /£ s = 0.118 ± 0.062 

and 

f£ u = 0.014 ± 0.003, f£ d = 0.036 ± 0.008, f£ s = 0.118 ± 0.062 , 



(4.5) 



(4.6) 
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FIG. 3: Dark matter- nucleon elastic scattering cross section for the parameter sets in Table HI 
(stars) and Table [Till (t riangles). The solid line is the current bound from CDMS Soudan 2004- 



2009 Ge 



26l |. The dashed line represents the projected bound from SuperCDMS Phase A. The 



dotted line represents the projected reach of the LUX LZ20T experiment, assuming 1 event sen- 
sitivity and 13 ton-kilodays. The graph is obtained using the DM Tools software available at 



http: / / dmtools.brown.edu 



while is defined by 



JTg 



1 



q=u,d,s 



(4.7) 



We can approximate f p ps /„ since fr s is larger than other /t 9 's and fr g - For the purpose 
of comparing the predicted cross section with existing bounds, we evaluate the cross section 
for scattering off a single nucleon, which can be approximated 



O,, 



m rfp 
7T 



(4.8) 



where m r is nucleon-dark matter reduced mass l/m r = l/m n + l/m^. Our results are 
shown in Fig. El for the parameter sets given in Tables [III and IIHI The predicted cross 

26l | for dark matter masses between 1 and 



sections are far below the current CDMS bounds 



4 TeV. However, there is hope that the model can be probed by the future LUX LZ20T 



experiment 



29]. 
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V. CONCLUSIONS 



We have presented a new TeV-scale model of decaying dark matter. The approximate 
stability of the dark matter candidate, ip, is a consequence of a global U(l) symmetry that 
is exact at the perturbative level, but is violated by instanton-induced interactions of a non- 
Abelian dark gauge group. The instanton-induced vertex couples the dark matter candidate 
to heavy, exotic states that mix with standard model leptons; the dark matter then decays 
to i + i~v final states, where the leptons can be of any generation desired. We have shown 
that a lifetime of ~ 10 26 s, which is desirable in decaying dark matter scenarios, can be 
obtained for perturbative values of the non-Abelian dark gauge coupling. In addition, by 
studying dark matter annihilations through the Higgs portal, we have provided examples 
of parameter regions in which the appropriate dark matter relic density may be obtained, 
assuming dark matter masses that are consistent with fits to the results from the PAMELA 
and Fermi-LAT experiments. The nucleon-dark matter cross section in our model is lower 
than the present bound from CDMS, but may be probed in future experiments. It might 
also be possible to probe the spectrum of our model at the LHC. 

The model in this paper provides a concrete, TeV-scale scenario in which dark matter 



decay is mediated 
dark gauge groups 



by 



instantons, and gives a new motivation for the study of non-Abelian 



30] . However, it is by no means the only possible model of this type. One 
might study variations of the model in which different annihilation channels are dominant, 
or the dark matter is lighter, or the standard model leptons are directly charged under 
the new non-Abelian gauge group. It may also be worthwhile to consider how low-scale 
leptogenesis and baryogenesis might be accommodated in this type of scenario. While we 
have assumed parameter choices motivated by the observed cosmic ray positron excess, one 
might incorporate the present model in a multi-component dark matter scenario if this were 
required to explain new results from ongoing and future direct detection experiments. 
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Appendix A: Mass mixing example 



In Sec. m we presented a diagrammatic representation of the mixing that takes the 
X states to standard model leptons. Here we study the numerical diagonalization of the 
corresponding fermion mass matrices, to demonstrate that mixing angles of the size assumed 
in our analysis are easily obtained. To simplify the discussion, we focus on mixing with 
standard model leptons of a single generation, which we denote by e and v. We include (1) 
Dirac masses for the \ fields: 

^Eh ^ {Hd)x « r + h ^ {Hd)x ^ r + q ^ {5d)x - r + di ^ ^ x< M + h - c - ' (A1) 

i 

where Hp = io 2 H* D . These terms generate a completely general two-by-two Dirac mass 
matrix for the % fermions. (2) Mixing between the \ fields and standard model leptons: 

C D gi(v)Xme C R + 92(v)Xu ) R e C R + K~L(H)e R 

+ 93(v)XdR u R + 9MXurVr + \vL(H)u R + h.c. 
(3) Mixing involving the vector-like leptons El and Er\ 

C D g 5 {v) xf R E L + g 6 (v) XurEl + M e E l Er + g 7 L(H)E R + h.c. 
We now write down the mass matrices which follow from Eqs. flAl|A2IIA3j) . For the neutral 
states, we work in the basis fl = (XuL XdL u r) and Jr = (x2?> v l)- The neutral mass 
terms can be written as /"M /}| + h.c, where 

/ A _„,_ n \ 



(A2) 



(A3) 



c 2 v D d 2 v D 

a 2 v D b 2 v D , (A4) 

^ 94V V g 3 v v y/2m u j 

assuming, for simplicity, that the vevs and couplings are real. Similarly, the mass terms 
for the charged states may be written f£M c f^ + h.c, where we assume the basis f£ = 
(XST, A xS, A El, e L ) and f~ R = x » x g>, E R , e R ). In this case, 



( c\Vd cliVd 
diVp biVD 

c 3 v D d 3 v D 

a 3 v D b 3 v D 

\ 



92Vr, 
9lVr, 

o 



gev v g5V v V2M E 
g 7 v \f2m e J 










(A5) 
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Given a choice of parameters, it is now a simple matter to compute the relevant mixing 
angles numerically. As an example, let us work in units of the dark scale vd, which we will 
assume is 4 TeV. In addition we take v v = vjj, Me = 1.5 vp and set the standard model 
lepton masses to zero (the conclusions do not change if we require realistic standard model 
lepton masses). If one assumes that only the following parameters are nonzero: 

{61, ci, b 2 , c 2 , 63, c 3 , g u g 2 , g 3 , 9a, 95, 9e, 97} = 

{1.9, 1.8, 1.8, 1.7, 2.1, 2.0, 0.02, 0.02, 0.02, 0.02, 0.7, 0.6, 1.0} , (A6) 

then one finds 

x£2 = 0-011 e Ro + ■■■ x2 = 0.011 e Ro + ■■■ 

xS = 0-012 v R0 + • • • X f L = 0.011 u% + ■ - ■ 

X { ul = 0.009 e L0 + ■ ■ ■ xtl = 0.010 e Lo + ■ ■ ■ 

where the fields on the right represent mass eigenstates. In addition, the non-zero mass 
eigenvalues are all larger than the ip mass if rrif < 1.2 vd, so that only decays to standard 
model leptons via the instanton vertex are kinematically allowed. Given the number of free 
parameters involved, one sees that the mixing angles are highly model dependent and can 
be easily set to the values assumed in Sec. [TTl 
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